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ABSTRACT 

The horizontal branch (HB) morphology of globular clusters (GCs) is most 
strongly influenced by metallicity. The second parameter phenomenon, first de- 
scribed in the 1960's, acknowledges that metallicity alone is not enough to de- 
scribe the HB morphology of all GCs. In particular, astronomers noticed that 
the outer Galactic halo contains GCs with redder HBs at a given metallicity than 
are found inside the Solar circle. Thus, at least a second parameter was required 
to characterize HB morphology. While the term 'second parameter' has since 
come to be used in a broader context, its identity with respect to the original 
problem has not been conclusively determined. Here we analyze the median color 
difference between the HB and the red giant branch (RGB), hereafter denoted 
A(V — J), measured from Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) photometry of 60 GCs within ~20 kpc of the Galactic Center. 
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Analysis of this homogeneous data set reveals that, after the influence of metal- 
licity has been removed from the data, the correlation between A(V — I) and age 
is stronger than that of any other parameter considered. Expanding the sample 
to include HST ACS and Wide Field Planetary Camera 2 (WFPC2) photometry 
of the 6 most distant Galactic GCs lends additional support to the correlation 
between A(V — I) and age. This result is robust with respect to the adopted 
metallicity scale and the method of age determination, but must bear the caveat 
that high quality, detailed abundance information is not available for a signif- 
icant fraction of the sample. Furthermore, when a subset of GCs with similar 
metallicities and ages are considered, a correlation between A(V — I) and cen- 
tral luminosity density is exposed. With respect to the existence of GCs with 
anomalously red HBs at a given metallicity, we conclude that age is the second 
parameter and central density is most likely the third. Important problems re- 
lated to HB morphology in GCs, notably multi-modal distributions and faint 
blue tails, remain to be explained. 

Subject headings: globular clusters: general 



Introduction 



It has been clear for decades that metallicity is the most influential factor governing the 
HB morphologies of Galactic GCs: metallicity is the fir st parameter. T h e earliest photo- 
graphic color- magnitude diagrams (CMDs) of GCs by, e.g. lArp et al.l (119521 ) ; ISandagd (119531 ). 
revealed that HB stars in metal-rich GCs tend to lie on the red side of the RR Lyrae insta- 
bility strip while HB stars in metal-poor GCs lie primarily on the blue side. As more an d 



more CMDs were assembled, exceptions to this rule were uncovered. ISandage fc Wallerstein 
( 119601 ) noted that M 13 and M 22 display HB morphologies appropriate for metal-poor GCs, 
despite the fact both GCs w ere of intermediate metal licity, and suggested that a difference 
in age might be responsible. ISandage fc Wildeyl (119671 ) noticed that the HB of NGC 7006 is 
redder than either M 13 or M 3, despite the fact that all three GCs appeared to have very 
similar metallicities. Such anomalies suggested the need for a second param eter that could 
account for differences in HB morphology not obviously caused by metallicity. Ivan den Bergh 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555, under program 
GO-10775. 
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( 119651 ) summarized the problem, stating that metallicity is not sufficient to explain the ex- 
tent of observed HB morphologies but that differences in age or He enrichment could explain 
the observed variatio ns. Both of these suggestions remain valid up to the present time. 
van den Berghl (119671 ) analyzed the integrated colors of 49 GCs with UBV photometry and 
concluded that 'at least two parameters (one of which is metal abundance) are required to 
describe globular clusters.' 

The second parameter p henomenon took on a greater significance with the seminal 
work of ISearle &: Zinnl (119781 ). Searle & Zinn recognized that GCs with unusually red HBs 
are relatively rare in the inner regions of the Galactic halo (Galactocentric radius, Rqc ^ 
8 kpc) but become increasingly common at greater Rqc (see their Figure 10). Searle & 
Zinn used this fact to argue that the inner halo was assembled early and in a fairly short 
time while the outer halo was assembled over an extended period of time. Current galaxy 
formation scenarios envision the outer regions of the Galaxy as the accumulated debris of 
the many accretion events that shaped the early evolution of the Galaxy. In that context, 
understanding the origin and existence of age and metallicity gradients in the Galactic GC 
population is as relevant now as it was at the time when Searle & Zinn first introduced their 
halo formation scenario. 

Subsequent efforts relating to the second parameter problem fall into two general cate- 
gories: those that attempt to measure the age difference between two (or a few) GCs with 
similar metallicities but markedly different HB morphologies and those that investigate the 
ensemble properties of a large sample of GCs. 

One canonical pair in the former category is NGC 288, with a blue HB, and NGC 362, 
with a red HB. The firs t CCD-based, differential photometric study of these clusters was con- 
ducted bv lBolte! Jl989h. Bolte conc l uded t hat NGC 288 is ~3 Gyr older t han NGC 362. The 
subseq uent work of Green fc Norris ( 199ol), Saraiedini fc Demarquel Jl990 ). and VandenBerg et al 
( 1990 ) echoed this result. By contrast, VandenBerg &: Durrell ( 1990T) reached a very differ- 
ent conclusion. Using the brightest RGB star in each GC 



VandenBerg fe Durrelll (I1990f ) 



corrected the photometry of NGC 288 and NGC 362 for their relative distances and found 
their main sequence turnoffs to b e roughly coincident, thereby suggesting a negligible age 
difference. iBellazzini et~aD (hoOlh found NGC 288 to be 2±1 G yr older than NGC 362 using 
three different techniques; in a follow-up, ICatelan et al.l (j2001al ) found that an age difference 
of 2 Gyr was plausible if both GCs have [Fe/H] ~ — 1.2 but that their synthetic HB models 
were unable to match the detailed HB morphology of either GC using canonical assumptions 
of average mass loss and dispersion on the HB. In considering differences between second 
paramete r pairs, it is important to understand just how similar the abundances are in both 
clusters. IShetrone &: Keand (120001 ) performed a detailed comparison of abundances in red 
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giants with 13 stars in NGC 288 and 12 in NGC 362. These authors concluded that NGC 288 
has a lower [Fe/H] than NGC 362 by 0.06 dex, that the average [a/Fe] ratios are very similar, 
and that both GCs exhibit variations among O, Na, and Al. 

The other canonical second parameter pair is M 13, with a blue HB, and M 3, with 
an intermediate HB. M 13 and M 3 are ~0.2 dex more metal poor than NGC 288 and 



NGC 362. IVandenBerg et al.l (119901 ) found evidence for an age difference but the quality 



of their CMDs did not allow a more definitive statement. Catelan fc de Freitas Pacheco 



( 119951 ) estimated an upper limit to the age difference of ~3 Gyr and concluded that it was 
insufficient to explain the H B morphologies, a ssumi ng both GCs have the same chemical 
composition. Alternatively, iJohnson &: Boltd (119981 ) suggested the difference in their HB 
morphologies was due to a difference in their He abundances. Johnso n fc Bolte suggested 
that M 13 had a He mass fraction ~0.05 greater than M 3. However, ISweigart fc Catelan 



( 119981 ) showed that a difference of AY ~ 0.05 at constant [Fe/H] would make the level of 
the HB brighter b y ^0.2 mag . Such a difference between M 3 and M 13 was ruled out by 
the photomet ry of iRev et aJj (120011) . who also concluded that M 13 is older than M 3 by 
1.7±0.7 Gyr. ISneden et al.l (120041 ) compared spectra of 28 red giants in M 3 with 35 in M 13 
and found the two G Cs' mean [Fe/H] values the same within their 1-a error bars. However, 
Sneden et al.l (120041 ) reported differences in light element abundance variations, particularly 
O, for which M 13 displays a larger range of variation than M 3 by ~0.5 dex. 

A discussion of the search for the secon d parameter in GC-to -GC comparisons would no t 
be complete without including the work of IStetson et al.l (119991 ) and iDotter et al.l (l2008bl ). 
These investigations used HST WFPC2 ph otometry to measur e the ages of the outer halo 
GCs Palomar 3, Palo mar 4, and Eridanus (IStetson et al.l 119991 ) and AM-1 and Palomar 14 
( Potter et al. 2008b ) relative to inner halo GCs with similar metallicities, M3 and M5^|. 
Each of the five outer halo GCs has a redder HB morphology than its comparison inner 
halo GC. Both studies concluded that the outer halo GCs are ~1.5-2 Gyr younger than M3 
and M5, provided the chemical compositio ns of the ou t er ha lo GCs are comparable to their 
inner halo counterparts. A recent study by lKoch et al.l (120091 ) found the abundances of Pal 3 
are essentially the same as inner halo GCs of similar metallicity and thus the relative age 
comparison with M 3 by IStetson et al.l (Il999l ) was justified. Overall, the outer halo GCs' 
che mical compositions r emain poorly u nderstood by comparison with the inner halo, see 
e.g. lPritzl et al.l (120051 ). ICatelanl (120001 ) found the reported age differences between Pal 4 
and Eridanus, with red HBs, and M 5, with an intermediate HB, too small to explain the 
difference in HB morphologies unless all three GCs are younger than 10 Gyr, assuming 



2 There is substanti al evi dence that M 5 is ac tually an outer halo GC currently near its perigalacticon, 
see Scholz et al. ( 1996 ) and Dinescu et al. ( 19991 ). 
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standard assumptions of mass loss on the RGB. ICatelan et al.l (l2001bl ) concluded that it 
was possible to explain the difference in HB morphology between Pal 3 and M 3 if the 
former has less HB mass dispersion and is younger than the latter. These examples indicate 
that if age is the second parameter, then it is our lac k of underst anding of mass l oss that 
confus es GC-to-GC comparisons as first pointed out by lRoodl (119731 ). The review by lCatelan 
( 120091 ) includes a thorough discussion of different mass loss prescriptions and their efficac 
Given the HB morphologies and relatively young ages of the most distant outer halo GC 
it is important to include them when considering the properties of the entire Galactic GC 
population, especially in the context of Searle & Zinn's halo formation scenario. 

The second category of second paramet er studies are t hose that consider the prop- 
erties of a large sample o f Gala ctic GCs, e.g. lSearle fc Zinnl (119781 ). In the first CCD-era 



study, ISaraiedini fc Kind ( 



989h compiled properties for 31 Galactic GCs. Among other 



things, ISarajedini fc King! (119891 ) examined the variation of HB morphology with age for 



GCs in a narrow range of metal abundance and showed that GCs with red HBs are signif- 
ic antly younger than tho s e with blue HBs by as muc h as 5 Gyr. The subse quent studies 
of IChaboyer et al.l (119921) . ISaraiedini et al.l (|1995T ). and lChaboyer et ajj ( 119961 ) updated and 
reaffirmed this result. Rosenberg et al. ( 19991 ) found that the GCs with Rqc > 8 kpc are, on 
average, younger than t hose with Rgr , < 8 kpc, which is consistent with age as the second 
parameter. By contrast, iRicher et al.l (119961 ) examined the CMDs of 36 GCs, found an age 
dispersion of ~1 Gyr with no significant age gradient in the Galactic halo, and concluded 
that this age range was too small for metallicity and age alone to explain HB morphology 
in GCs. 

Theoretical efforts have provided further insights into the complexities of HB mor- 
pholo g y. In p articular, s evera l studies have applied the synthetic H B mo del develo p ed by 
Rood! Il973h . iLee et all Jl99ph . ICatelan fc de Freitas Pachecol Jl993h . and iLee et all Jl994T ) 
used synthetic HB models to explore the interplay of HB morphology, metallicity, and age 
in the CMD. These studies demonstrate d considerable degeneracy in t h e HB morphology- 
metallicity diagram and, in particular, ICatelan &: de Freitas Pachecol (119931 ) argued that 
unless absolute values of the He abundance, [a/Fe] ratios, and RGB mass loss were known , 
HB morphology did not constitute a reliable age constraint. Nevertheless, ILee et al.l (119941 ) 
concluded that there was evidence for an age dispersion of < 5 Gyr among Galactic GCs of 
similar metallicity but markedly different HB morphologies. A number of theoretical studies 
have also focused on second parameter pairs or triads but, unfortunately, the lack of a firm 



3 The exception in the outer hal o is NGC 2419. A num ber of dedicated HST ph otometric studies have 
focused on this massive, distant GC. lHarris et alj (|l997l ) and lSandquist fc Hessl (|2008h both found NGC 2419 
to be an old, metal poor GC with a blue HB. 
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theoretical und erstanding o f mas s loss along the RGB has impaired the these efforts, as first 
pointed out by [Rood Jl973h . See ICatelanl (120091 ) for a recent review. 



As it pertains to HB morphology, the use of the term 'second parameter' has expanded 
to cover a broad range of comple x behaviors. In p articu lar, this includes the faint extent 
of the blue HB tail in some GCs. iFusi Pecci et al.l (119931 ) analyzed 53 GCs and found that 
the lengt h of the HB in the CMD -and the extent of the blue tail-is correlated with central 
density. iBuonanno et al.l (119971 ) also reported a link between an extended blue t ail and 



centra l density in GCs and concluded that 'environment is "a" second parameter'. iSmith 
(120041 ) showed that there is a correlation between central density and HB morphology for 
intermediate met allicity GCs (—1.7 < [Fe/H ] < —1.3) where the second parameter effect is 



most pronounced. iRecio-Blanco et al.l (120061 ) analyzed the properties of 54 Galactic GCs with 



homogeneous photometry from an HST WFPC2 Snapshot Survey ( Piotto et al. 20021) and 



concluded that more massive GCs tend to have more extended blue HBs. IRecio-Blanco et al. 



( 120061 ) recognized a link between the effective temperature of the hottest HB star in a GC 
and its mass (as inferred from its integrated luminosity) and suggested that self-pollution 
could explain the existence of faint blue tails in preferentially massive GCs. 



E vidence abounds for th e presence of chemical abundance variations in all GCs (IGratton et al. 

2004bl ; ICarretta et al.lfiooaal Jbh and the possibl e correlation between the degree of abundance 
variations and the faint extent of the blue HB (ICarretta et al.ll20071 ). However, it is unlikely 
that the existence of faint, blue tails in the HBs of some GCs is directly related to the appear- 
ance of anomalously red HBs in others. It is also unlikely that chemical abundance variations 
will unduly affe ct age estimates of most GCs, provided the total metal content is constant 
across all stars (jPietrinferni et al.l 120091 ) . For cases in which the total metal abundance is 
not constant, or there are distinct stellar populations present in the CMD, age estimates 
are necessari ly more complicated and careful analysis of each of these GCs is needed. See 
Piottol (120091 ) for a recent summary of GCs with multiple stellar populations visible in the 
CMD, several of which were discovered with photometry from the ACS Survey of Galactic 
GCs. The complex issue of multiple stellar populations in GCs remains poorly understood 
and the extent to which multiple-population GCs permeate the Galactic GC population is 
unknown at present. 

To summarize, the presence of metal poor GCs with red HBs that primarily reside in 
the outer Galactic halo is well-known observationally. Although the present study is focused 
on the Galactic GCs, including some of those associated with the Sagittarius dwarf, there 
is ample eviden ce to suggest that metal poor GCs wit h red HBs are also found in the Mag - 
ellanic Clouds (|Johnson et al.l Il999l : iGlatt et al.l 120081 ) and Fornax (IBuonanno et al.l Il999l ) . 
Mackey fc Gilmord (120041 ) summarize our current knowledge of the GC populations in the 
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Galaxy and its satellites in the HB morphology-metallicity diagram. Despite an abundance 
of observational evidence, no consensus has been reached as to what parameter(s) are respon- 
sible for the appearance of relatively red HBs in metal poor GCs. Age is frequently offered 
as a candidate but considerable doubt still remains because the age difference claimed-or 
required by theoretical studies of HB morphology-is too large to satisfy the observations. 

The existence of a homog eneous database of de ep, high quality photometry from the 



ACS Survey of Galactic GCs (jSarajedini et al.l 120071 . Paper I in this series) has motivated 



a re-examination of HB morphology and its relation to a variety of GC properties. The 
paper is organized as follows: $2] describes the data sample; §3] describes the methods that 
were employed to determine the HB morphologies; §U describes the sources of GC ages 
and provides some discussion of complicating factors in age determination; §5] presents the 
analysis performed on the assembled data and discusses the important results; and, finally, 
£0 provides a summary of the salient points. 



The observational data 



2.1. The ACS Survey of Galactic GCs 



2008 



The photometric catalog of the ACS Survey of Galactic GCs (lAnderson et al 
Paper V) consists of 65 GCs observed in F606W and F814W. Most details concerning 
the data reduction and calibration are provided in Paper V but, since its publication, the 
data have bee n adjusted to account for updated HST/ACS Wide Field Channel zero-points 
(IBohlinl 120071 ); the n ew -F606iy and F814 W zeropoints are, respectively, 22 and 35 mmag 
fainter than given by ISirianni et al.l (120051 ) . 



Not all GCs in the ACS Survey catalog will be considered in the analysis performed 
in the following sections, u Cen (NGC 5139) and M 54 (NGC 6715) have been excluded 
because their complex CMDs and abund ances indicate that these objects are not GCs in 
the traditional sense. ISiegel et al.l (120071 ) analyzed the M 54 CMD and measured ages for 
the stellar populations present and also estimated the metallicities of the stellar populations 
for which no spectroscopic information was available. Because M 54 lies at the center of 
the Sagittarius dwarf galaxy, it is not a simple process to distinguish between the stellar 
population(s) that belong to M 54 and those of the Sagittarius dwarf galaxy, uo Cen has long 
been known to contain a ran g e of different metallici ties as well as multiple stellar populations 
in the CMD (ILee et al.lll999l ). IJohnson et al.l (120081 ) presented a spectroscopic analysis of 180 
red giants in u Cen that revealed at least four distinct metallicities and a total range of ~ 1.5 
dex in [Fe/H]. Other GCs for which multiple stellar population evidence exists will be noted 
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in the following sections. iPiottol (120091 ) provided a summary of the GCs known to exhibit 
multiple stellar populations in their CMDs circa 2008. Three more GCs were excluded: Pal 1 
and E 3 for lack of identifiable HB stars and Pal 2 due to extreme differential reddening (see 
Paper I). In total, 60 GCs from the ACS Survey are considered in the following analysis. 



2.2. The six most distant GCs 



The ACS Survey provides coverage of the Galactic GC population out to ~20 kpc. 
In order to give proper consideration to the outer Galactic halo, where the second param- 
eter phenomenon is most pronounced, the sample was exten ded to include the following 
GCs: NGC 2419 Jliarris et all Il997l: ISandauist fc Hessll2008h: P al 3, Pal 4, and Eridanus 
( Stetson et al. 19991 ): and AM-1 and Pal 14 ( Potter et al. 2008b ). These are the six most 
distant GCs in the Galaxy with 70 < Rgc ^ 120 kpc (jHarrislll996l ). Unfortunately, no pho- 
tometry of comparable quality is currently available for the 10 or so GCs that lie between 
these two groups (25 < Rqc ^5 50 kpc)B In order to distinguish between the homogeneous 
ACS Survey sample and the larger, heterogeneous sample, the statistical analysis in section 
$5] will be presented with and without the six outer halo GCs. 



2.3. Additional GC parameters 



In addi tion to HB morphologies ( $!]) an d ag es (Bj), several GC parameters were ex- 
tracted from lMackev &: yan den Berghl (120051 ). The lMackey fc van den Berghl ( 20051 ) catalog 
is based on the lHarrisI Jl996h catalog but includes several updated quantities and, impor- 



tantly, provides quantities in absolute physical units. These u n its ar e more appropriate for 
comparisons than the observational units given in the lHarrisI (119961 ) catalog. The quanti- 
ties obtained from the iMackey fc van den Berghl (120051 ) catalog are: Galactocentric distance 
(Rgc), in tegrate d abso lute V magnitude (My), half-light radius (Rh), and tidal radius (Rt). 
From the lHarrisI (119961 ) catalog we extracted central luminosity density (po)- Neither catalog 
includes basic measurements for Lynga 7. Inclusion of these quantities in the following anal- 
ysis makes it possible to explore the possibility that a GCs location in the Galaxy, structural 
properties, and total mass (assuming the mass-to-light ratio is constant) can influence its 
HB morphology. 



4 Pal 2 is one of these, and part of the ACS Survey catalog, but its CMD is so obscured by differential 
reddening that no useful information can be extracted from it (see Paper I). We note that during HST Cycle 
17, GO 11586 will use ACS to image three GCs at R G c ~ 40 kpc: Pal 15, NGC 7006, and Pyxis. 
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3. Measuring HB morphology 



In this study, HB morphology is characterized by the difference between the median color 
of the HB stars and the median color of the RGB at the level o f the HB, denoted A(V — I). 
This metric is essentially the same as d^-v fjSarajedinil Il999l ) and is less likely to become 
saturated th an t h e mos t frequently used HB morphology parameter, (B— R)/(B+V+R) in- 
troduced by iLed (119891 ) . As a corollary, the magnitude level of the HB in F606W has also 
been measured; it provides a reference for GC distance estimates as discussed below. 



3.1. The mean magnitude level of the HB 

The mean level of the horizontal part of the HB was measured from the ACS Survey 
data relative to the HB of M 5 (NGC 5904). M 5 has a well-populated HB that extends 
approximately 0.7 mag in F606W — F814W from the red side to the blue and, therefore, can 
overlay the HB of almost any other GC. Measurements were made by matching the CMD of 
each GC to that of M 5 in the region of the HB by making vertical and horizontal adjustments 
to the comparison CMD until it most closely overlaid that of M 5. In the majority of cases, 
the comparison GC had enough overlap to make the best alignment unambiguous. A few 
clusters with purely blue HBs, such as NGC 6254, present some difficulty because even the 
reddest HB stars do not become 'horizontal'. The measured HB levels are listed in Table 
CD Typical measurement uncertainties are ~0.05 mag and are dominated by the uncertainty 
inherent in aligning the HBs. The formal photometric errors of most HB stars are a few 
thousandths of a magnitude. The uncertainty is larger in cases where differential reddening 
is significant or the HB is otherwise difficult to quantify. As an example, Figure [1] shows HB 
levels for NGC 1261 and M 92 (NGC 6341). 



3.2. The color difference, A(V - I) 

Measurements were obtained by drawing an outline around the HB and the RGB at 
the level of the HB and then calculating the statistics of the stars within each outline. This 
process is illustrated in Figure [1] for two GCs, NGC 1261 and M 92. RR Lyrae stars are 
not given special consideration in this analysis: they are included if they are present in the 
ACS observations. For the ACS Survey data, all determinations of HB membership were 
performed in the ACS F606W — F81AW CMD as shown in Figured! In order to allow direct 
comparison with GCs observed in different bands, the color difference presented in Table [1] 
is given in terms of V — I rather than ^6061^ — F81AW. This was achieved by transforming 



-li- 



the ACS magnitudes of each individual HB star into V and I using the equations provided 
by lSirianni et al.l (120051 ) and then calculating the median V — I values of the HB and RGB. 
The median was chosen to represent the central tendency of A(V — I) because the HB stars 
are not normally distributed in color space. Hence, the more commonly-used mean and 
standard deviation are not the most appropriate choices. If the HB stars were normally 
distributed in color space, then the mean and median would be equal. In fact, about one 
half (one quarter) of the GCs in the sample have means and medians that differ by more 
than 5% (10%). For GCs with intermediate HBs, such as M 3 and NGC 1261, the mean and 
median are substantially different. 

Table [1] provides the quantities used to determine A(V — I)- In addition to the median 
V — I of the HB and RGB, 1-er errors on these quantities are listed and the quadrature sum 
is given as the uncertainty on A(V — I). Uncertainty in the median V — I was estimated 
by bootstrapping with replacements performed 10,000 times on both the HB and RGB of 
each GC. The error bars on the median represent the range of V — I within which 68% 
of bootstrapped medians lie. To compliment the median and uncertainty, we also give 
the mean absolute deviation (MAD) for the HB and RGB. The mean absolute deviation 
represents the intrinsic spread in each individual data set. Consider, for example, the HBs 
of 47 Tuc (NGC 104) and M 3 (NGC 5272). Both have well-populated HBs and thus small 
uncertainties, o < 0.01 determined from bootstrapping, but 47 Tuc has a tightly clustered 
red clump of HB stars and a small MAD=0.023 while M 3 has a more broadly distributed 
HB and a larger MAD=0.208. 



3.3. Comparison with other HB morphology parameters 



Perhaps the most frequently used HB morphology metric in the literature is (B— R) / (B+V- 
the difference between the number of blue HB stars (B) and the number of red (R) normal- 
ized by the total number of HB stars including variables (V). Its major inconvenience is that 
it saturates if all stars lie on one side or th e other of the instability strip. T he left panel of 
Figure [2] compares (B— R)/(B+V+R) from iMackev fc van den Berghl (120051 ) and A(V — I) 
from this paper for the 66 GCs in the present sample. There is an obvious correlation 
between the two but A(V" — I) continues to vary after (B— R)/(B+V+R) saturates. 



The most recent large-scale study of HB morphology, that of iRecio-Blanco et al.l (120061 ). 
used the maximum effective temperature encountered along the HB [Log T R ff(HB)1 deter- 



mined with the use of theoretical zero-age HB sequences from ICassisi et al.l (119991 ) . The 



Log T e fj(HB) metric appears to be complimentary to A(V — I) or (B— R)/(B+V+R) be- 
cause, while it lacks sensitivity through the middle, it is more sensitive at the extremes. This 



R): 
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is demonstrated in the right panel of Figure [2j In the absence of an ideal HB morphology 
metric that quantifies every feature of interest, it is beneficial to identify a metric that is 
well suited to a particular problem. 



4. Globular cluster ages 



4.1. Relative ages 



Marin-Franch et al.l (120091 . Paper VII) presented relative ages measured from the ACS 



Survey data. These ages were derived by first measuring the absolute magnitudes of the 
main sequence turnoffs (MSTOs) for all GCs and then interpolating in isochrone-based grids 
of MSTO as a function of age and metallicity. The same analysis was performed with four 
different isochrone libraries and the distribution of relative ages with metallicity was shown 
to be independent of the isochr o ne lib rary. For the purpose of the pr esent paper, r e lative 



ages derived from iDotter et al.l (120071 . Paper II) isochrones using the IZinn fc Westl (11984 



hereafter ZW84) metallicity scale have been adopted. Furthermore, we assumed that GCs 
with [Fe/H] < -1 have [a/Fe]=+0.3 and GCs with [Fe/H] > -1 have [a/Fe]=+0.2. The 
ages measured in this fashion and with these assumptions were placed on a relative scale by 
dividing out the average age of the most metal-poor GCs ([Fe/H] < —1.8), in this case 13.3 
Gyr. In the present study, this factor has been retained, resulting in absolute ages for the 
particular isochrone library, metallicity scale, and method of age determination employed in 
Paper VII. 



4.2. Isochrone fitting 



In order to put the ACS Survey ages on the same scale as the six outer halo GCs, 
and to demonstrate that later results do not depend on the method used, age estimates 
were determined using isochrone fitting to the CMDs of all GCs in the data set. The fit s 



were pe rformed using isochrones from Paper II for the ACS Survey data and IDotter et al. 



(j2008al ) for the outer halo GCs. The same zeropoint corrections applied to the ACS data, 
as described in §2.1[ were applied to the isochrones in the ACS photometric system as well. 
The underlying luminosities, t emperatures, s urface gravities, and color transformations of 
the isochrones in Paper II and IDotter et al.l (j2008al ) are identical. The only difference be- 
tween the two is the photometric system. The outer halo GCs were measured in F555W 
and F814W, whereas the ACS Survey data were taken in F606W and F814W. Several 
GCs were excluded from isochrone fitting because they are known to harbor multiple stellar 



13 



populations. The presence of more than one stellar population makes a single age deter- 
mination insuffici ent to characterize the GC. Ages were n ot determined for the following 
GCs: NGC 1851 iMilone et alihoosl Paper III), NGC 2808 JPiotto et alil2007h . NGC 6388, 
NGC 6441@, and NGC 6656 <jPiottoll2009h . 



It is likely that other GCs with multiple populations exist, even within the ACS Survey 
data, but to have avoided detection thus far the separation between the populations must 
either be so small that they overlap to a great extent or that reddening has obscured the 
separation in the CMD. A second population may also represent a small fraction of stars and 
have little or no apparent influence on the HB morphology The recent discovery of a second 
subgiant branch in 47 Tuc (NGC 104) by lAnderson et al.l (120091 ) is an example: this second 
population accounts for only ~10% of the stars in the core. Further out, th e CMD was too 



sparse ly populated to detect this population. It should also be noted that lAnderson et al. 



( 120091 ) found the main population of 47 Tuc to have a spread in the CMD greater than 
can be accounted for by the photometric errors. The authors fou nd it unlikely that this 
spread was due to binaries, differential reddening, or a depth effect. lAnderson et al.l (120091 ) 
explored the possibility that the color spread could be caused by an intrinsic varia tion in the 
He abundance or [Fe/H]. Another case is M 4 (NGC 6121). iMarino et al.l (120081 ) presented 
spectroscopic evidence for a bimodal abundance distribution along the RGB of M 4 as well 
as a broadening of the RGB in U — B. Whether M 4 contains two distinct stellar populations 
or a continuous distribution is not clear at present. If it is the former, the difference is small 
en ough to remain und etected on the main sequence in the exquisite photometry presented 
by iBedin et al.l (120091 ) . If it is the latter, age determinations are expected to be reasonable 
(see the discussion at the end of §4.31) . 

The procedure used to determine ages involved a two step process. In itial est i mates 
for [Fe/H], distance modulus, and reddening estimates were taken from the iHarrid (119961 ) 
catalog, 2003 revision. Estimates of [a/Fe] were chosen as appropriate for a given value 
of [Fe/H]. All efforts were made to follow these values as closely as possible with minor 
adjustments made to improve the fit to the unevolved main sequence first and the RGB 
second. In some instances, isochrones at the catalog values were not able to provide an 
adequate fit to the CMDs and, in these cases, the initial estimates of [Fe/H], distance, 
and reddening were allowed to vary until the fit obtained was acceptable. This step was a 
necessity in a number of cases because the published values for some lesser-known GCs are 
uncertain. 



5 Given the great s i milari ties between the CMDs of NGC 6388 and NGC 6441 we have chosen to exclude 
both although Piotto (2(309) only showed evidence for multiple populations in NGC 6388. 
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For example, iDotter et al.l (l2008bl ) found that the main sequence and RGB morphologies 
of AM-1 are comparable to those of M 3 although the former is listed at [Fe/H] = —1.8 
while the latter is listed at [Fe/H] = —1.57 in the Harris catalog. Another example is the 
distance modulus of NGC 6254: the Harris catalog value is DMy = 14.08 but the present 
study produces a value that is ~0.2 mag fainter. This result is confirmed by both isochrone 
fitting to the main sequence and the level of the HB from Table [TJ These results are not 
meant to diminish the value of the GC catalogs, which are invaluable resources, but merely 
to emphasize that there are still many GCs in the Galaxy whose properties are poorly 
constrained. The models, which are homogeneous in terms of physical assumptions and 
ingredients, can prove useful in determining the relative differences between two GCs if one 
has well- measured properties and another does not. 

Several of the GCs in the ACS Survey are heavily obscured by interstellar reddening. 
Differential reddening confuses the age determinations in such GCs because the reddening 
line is nearly perpendicular to the most age-sensitive features: the main sequence turnoff and 
the subgiant branch. Differential reddening causes these features to be spread out and, in 
light of this trend, we have attempted to correct for its effects in the following manner. The 
first step involves the construction of a fiducial sequence for the GC. The fiducial sequence 
is centered approximately at the MSTO. More stars lie below the MSTO than above but 
the differential corrections are weighted based on the angle between the stellar sequence 
and the reddening line. Since the reddening line is almost perpendicular to the subgiant 
branch, the differential correction is approximately equally weighted by stars above and 
below the MSTO. Care is taken to avoid unresolved binaries. From the fiducial sequence, 
each individual star yields a color residual, taken along the reddening direction in the CMD. 
From these residuals a reddening map is created by finding the median residual in each 
256 x 256 pixel square of the image. Then each star is corrected, along its reddening line, 
by an amount that is interpolated from the 16 x 16 points of the reddening map. A full 
explanation of these procedures, the differentially corrected CMDs, and reddening maps will 
be provided in a forthcoming paper (I. King, 2009, in preparation). 

The best-fit age of a given GC was estimated by determining the isochrone that best fit 
the CMD from the MSTO through the subgiant branch. The uncertainty was derived from 
the intrinsic scatter in the CMD and/or the inherent mismatch between the models and the 
data. The range of ages that allowed the isochrones to envelope the bulk of the stars at the 
MSTO and on the subgiant branch are taken to be l-er uncertainties. The age uncertainty is 
only based on the fitting procedure described and does not account for uncertainties in the 
input physics of the stellar models or differences between the chemical composition assumed 
in the models and those actually present in the stars (for more on this point, see §4.3jl . 
Hence the uncertainties include the random component but exclude the systematic. Still, it 
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is worthwhile to consider that incomplete knowledge of chemical composition and incomplete 
treatment of the physics, such as rotation and convection, in the current generation of 1-D 
stellar evolution models are primary sources for the systematic error present in the analysis. 

The results of isochrone fitting are listed in Table [2j If a differential reddening corrected 
CMD was employed in the fit, an asterisk (*) appears after the name in Table [2j Figures [3] 
through [5] demonstrate how the isochrone fits were achieved for three different GCs. Figure 
E] shows a case where the isochrones trace the stellar population throughout the CMD (left 
panel) and how the age, 12.5±0.5 Gyr, was measured (right panel). Figure @] shows the 
differential reddening-corrected CMD of NGC 3201. The corrected CMD reduces the scatter 
about the subgiant branch but the age uncertainty is still larger than in the case of NGC 6362. 
Figure shows the worst case scenario where the reddening is low and the isochrones match 
the unevolved main sequence and RGB but not the shape of the age sensitive region. The 
slope of the models is shallower than the data along the subgiant branch and, although the 
data define a narrow sequence, the age is poorly constrained. This case is likely to be severely 
affected by the systematic errors described in the preceding paragraph. 

As a consistency check, it is useful to plot the absolute magnitude of the HB as a function 
of metallicity. The absolute magnitude of the HB is obtained by subtracting the F606W 
distance modulus derived from the isochrone fits (Table [2]) from the apparent magnitude of 
the HB (Tabled]). The resulting quantity is plotted in Figure [6] for the ACS Survey clusters. 
The points define a relatively tight relationship. A linear, least squares fit to these data for 
GCs with [Fe/H] < -1 gives: 

<M F6mw (HB) >= (0.227 ±0.011) [Fe/H] + 0.802 ±0.020, (1) 

which is shown as the solid line in Figure EJ The quoted errors in equation (jTJ are due only 
to scatter in the data and do not include measurement uncertainties. The fit is limited to 
clusters with [Fe/H] < —1 so as to exclude clusters with predominantly red HB morphologies. 
The justification for restricting the fit is that red HB GCs often lack RR Lyrae variable stars. 
Some GCs with completely blue HBs were included in the fit even though they may not have 
RR Lyrae stars. As described in §3.1[ the uncertainty in the apparent magnitude of the HB 
is approximately 0.05 mag per GC and this is shown by the two dashed lines. All but a 
handful of points lie within the dashed lines. 

If equation ([T]) is transformed from F606W to V using the synthetic color transforma- 
tions employed in Paper II, the slope increases by < 0.01 because of a very slight metal- 
licity dependence of V — F606W on [Fe/H] and, assuming a character istic temp e rature 



for RR Lyrae stars of LogT=3.83, the intercept increases by 0.09 mag. IChaboyerl (Il999l ) 



gave M V (RR) = (0.23±0.04)([Fe/H]±1.6) ± (0.56±0.12). Other recent estimates include 
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M v (g.B)=(0.22±0.05)([Fe/H]+1.5)+(0.56± 1 07) flGratton et al. 120031) and a further refine- 
ment of the slope to 0.214±0.047 mag/dex ( Gratton et all 120041 ) . If equation ([!]) is trans- 



formed again so that the independent variable is ([Fe/H] + 1.6) and the dependent variable 
is M v rather than M F606W , it becomes M V (HB) = 0.235([Fe/H] + 1.6) + 0.53 where the 
error bars have been omitted for brevity. The relationship between My(HB) and [Fe/H] 
derived here is within the l-er uncertainties from each of the three determinations previously 
mentioned. 



4.3. Metallicity scales and chemical abundance variations 



Figure [7J compares [Fe/H] from Table [2] and the ZW84 scale with the iKraft fc Ivans 



(120031 . hereafter KI03) scale for 47 GCs in common between KI03 and this paper. (KI03 
provided [Fe/H] values measured with three different model atmospheres; this paper adopts 
their measurements based on MARCS models.) For reference, the line of equality is drawn 
on bo th panels of Figure [7J Paper VII compared ages on the ZW84 and ICarretta fc Gratton 
( If997l . hereafter CG97) scales. As discussed by KI03, the CG97 scale is consistently higher 
than either the ZW84 or KI03 scales, by 0.2-0.3 dex in general. It is beyond the scope of 
this paper to determine which scale is superior to the others^]. Suffice it to say that ages 
and metallicities from Table [2] and Paper VII on the ZW84 scale will be used in the analysis 
that follows. However, speaking hypothetically, if it is assumed that the KI03 scale is the 
'correct' [Fe/H] scale, then the residuals from the solid line in each panel of Figure [7] can be 
interpreted as uncertainties in the adopted scales. We shall return to this hypothesis in §5.11 



The ages reported in Tables H] are based on isochrones that assume an a-enhanced 
composition. They do not explicitly account for variations in the He or C+N+O abun- 
dances but, if such variations are present, th ey should certainly influence age determinations 
dVentura et aD 120091 : b'Antona et all l2009h . While the extent of He and C+N+O vari- 
ations in GCs are largely unknown at present, the existence of light element abundance 
var iations is quite clear. T he review by iGratton et all (l2004bl ). the recent results presented 
by ICarretta et aD (ho09al lbl). and many other studies have shown the widespread presence 
of chemical abundance variations within individual GCs. These abundance variations are 
among the light elements (C, N, O, Na, Mg, and Al) and do not extend to heavier elements 
such as Fe in the majority of GCs. Models for GC self-enrichment predict that varying 



Recently, ICarretta et al.l (|2009cl ) presented a new GC metallicity scale based on the largest spectroscopic 



survey of GC red giants to date. The analysis reveals a very close agreement betw een the new me t allicity 
scale and that of KI03. This is in spite of the fact that KI03 used Fell lines while ICarretta et al. I (l2009ch 



relied on Fel lines. On the other hand, the new scale is consistently ~0.2 dex lower than that of CG97. 
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degrees of He variation sho uld accompany the other abundance variations (ID'Antona et al. 
20021 ; iDecressin et all 120071 ). The influence of such abundanc e variations on GC age s is not 
fully understood at present. However, the recent results from ICarretta et al. r|2009al Jbh lend 
some insight into the problem of deriving GC ages in the presence of light element abun- 
dance variations. Those authors identify three groups of stars within GCs based on their 
chemistry. They are: the 'primordial' stars or stars which have halo-like abundances, char- 
acterized by super-solar O and sub-solar Na; the 'intermediate' stars, with slightly less O 
and mor e Na than the p rimord ial stars; and the 'extreme' stars, characterized by [O/Na] 
< —0.9. Carretta et al. (|2009al Jbh state that all GCs (measured so far) contain primordial 
and intermediate stars but not all GCs contain the extreme component. Furthermore, the 
majority of GC stars are in the intermediate group. 



How do these findings influence GC age est imates? IPietrinferni et al.l (120091) compared 
stellar evolution models from the BaSTI library ( Pietrinferni et al. 20041. 2006T) with chemica l 
compositions representative of the primordial and extreme groups. IPietrinferni et all (120091 ) 
state that, as long as the total amount of C+N+O remains constant, it is safe to use a- 
enhanced isochrones to derive GC ages. However, if the total amount of C+N+O varies 
among the stars in a given GC, then it is necessary to use models with the appropriate 
CNO abundances to derive GC ages. This is all without consideration of a potential change 
in the He content that is predicted to accompany the light element abundance variations. 
The He abundance is important because models of chemical enrichment in GCs predict 
that the abundance variations from the primordial values should b e accompanied by an 
increase in He, whether the sourc e of the enrichment is massive stars (IDecressin et al.ll20071 ) 
or intermediate mass AGB stars (ID'Antona et al.l 120021 ) . The extreme stars should be most 
He-enriched, the intermedi ate stars shou l d be s lightly He-enriched, and the primordial stars 
should have primordial He. ISalaris et al.l (120041 ) investigated the He content of 57 GCs using 
the R parameter: the ratio of the number of HB stars to the number of RGB stars brighter 
than the HB. Within the uncertainties, they found little or no evidence for a spread in the 
He content o f GCs with (B-R)/(B+V+R) < 0.8. For GCs with (B-R)/(B+V+R) > 0.8, 
Salaris et al.l (120041 ) found a larger spread and a tendency toward higher He abundances. 
These authors state that the apparent trend towards higher He abundances in GCs with 
blue HBs may be related to increased evolutionary timescales for the lowest mass HB stars 
not entirely accounted for in their calibration of the R parameter or the genuine presence 
of He-rich stars in the blue HBs of some GCs. As long as the total metal content remains 
constant, a modest spread (AY < 0.05) in the He content of a GC will not sig nificantly alter 
the level of the main sequence turnoff or subgiant branch (IDotter et aJj 120091 ) and therefore 
not confuse the age determination. If the total metal content varies or the spread in He 
is larger, age estimates will require great care and detailed abundance information. It is 
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important to keep in mind the complexities of measuring GC ages in the context of chemical 
abundance variations-and their proposed origins-in GCs. 



4.4. Age-Metallicity Relations 

The average difference between the ages derived in §4.21 and in Paper VII is —0.104 Gyr; 
the standard error of the mean is 0.105 Gyr and the standard deviation is 0.781 Gyr. Figure 
[8] shows the normalized age differences of 56 GCs. The age differences are calculated by 
subtracting the Paper VII age from the age listed in Table [2] and dividing by the quadrature 
sum of the age uncertainties from both sources. In accordance with Paper VII, metallicity 



is represented 


in Fif 


(Salaris et al. 


1993) 



1-cr. The primary reason for disagreement in ages is most likely the adopted metallicity scale 
since the largest systematic deviation occurs around [M/H]~ —1 and the metallicity scales 
deviate there by as much as 0.3 dex (see Figure [7]). 

Figure M comp ares the AMRs from this paper with those of Paper VII (top panel) and 



VandenBergi (120001 ) (bottom panel). Despite the age differences already addressed in Figure 
El and the preceding paragraph, both AMRs reveal the same general features. The main 
difference is that the ages from Table [2] show smaller dispersion at low and high metallicities. 
The appearance of a separate trend beginning at [M/H]=— 1.5 and extending to Pal 12 and 
Ter 7 (the two youngest GCs) is nearly identical in both AMRs. The open circles in the 
upper panel are the outer halo GCs not present in the ACS Survey data; the se additional 



GCs s trengthen the trend already present in the ACS Survey data. Given that IVandenBerg 



( 120001 ) carried out his analysis using a heterogeneous collection of data in B — V and V — I, 
and used distances set by the level of the HB, it is encouraging that all three AMRs reveal 
the same basic trends. 



5. Results & Analysis 

5.1. A(V — I) and the second parameter 

HB morphology is most strongly influenced by the first parameter, metallicity. The 
metallicity- A (V — I) diagram is shown in Figure [TU1 For the ACS Survey data, the inner 
halo GCs (Rgc < 8 kpc) are plotted as circles and the outer halo GCs (Rgc > 8 kpc) as 
squares. The six most distant GCs are plotted as triangles and only appear in the left panel 
because they were not considered in Paper VII. The left panel shows the metallicity scale 
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from Table [2] and the right panel shows the ZW84 metallicity scale from Paper VII. The 
error bars shown are only due to measurement errors in A(V — I) as listed in Table [TJ 

In order to more clearly identify the second parameter, an attempt was made to remove 
the influence of metallicity from the HB morphology data, thereby exposing the second 
parameter effect. To accomplish this, a function was fit to the inner halo GCs which cover the 
full range of metallicity and exhibit relatively little scatter in the [M/H]-A(V — I) diagram, 
wit h the notable ex ception of NGC 6584. NGC 6584 is presently located at Rqc ~ 7 kpc 
but iLee et all (119941 ) noted its position in the metallicity-HB morphology diagram and large, 
positive helioce ntric velocity and sug gested it is actually an outer halo GC. This suggestion 
is supported by lDinescu et all (119991 ) who calculated its apogalacticon at 12.6 ± 2.4 kpc. As 
a result, M 5 (see and NGC 6584 were excluded from the fit. The fitting function is 
made up of two parts with x=[M/H], 



A(V-I) fit (x)=f(x)+g(x). 



(2) 



The first part, /, resembles a Fermi- Dirac fu n ction that transitions from blue to red as 
metallicity increases, as suggested by ICatelanl (120091 ) . The second part, g, is a quadratic 
that allows for the remaining curvature, in particular that the most metal-poor GCs turn 
back to the red with decreasing metallicity as can be seen in Figure [101 The two parts are 



f(x) = a - ai 



exp 



and 



g(x) = b + h\ x + b 2 x 2 



(3) 



(4) 



The a's and 6's were determined from a least squares fit. For comparison, and to demonstrate 
the robustness of the method, the same functional form was fit to the inner halo A(V — /) 
data using the metallicity scale from Table [2] (shown in the left panel of Figure [TOT) and the 
ZW84 metallicity scale from Paper VII (right panel). The fitting function coefficients are 
reproduced in Table [31 

The search for correlations focuses on the difference between the measured A(V — I) 
of an outer halo GC and the A(V — I) value of the inner halo trend (as represented by the 
fitting function) at that GCs metallicity. This quantity will henceforth be referred to as 
A (Fit — Measured). Figure [TT1 shows how A (Fit — Measured) varies with metallicity; the 
symbols are the same as in Figure [101 With the exception of NGC 6584 the inner halo 
GCs exhibit residuals of ~0.1 or less. Altogether, metal-rich GCs ([M/H]> —0.8) show very 
little scatter in Figure [TTJ Evidently, metallicity alone is almost sufficient to describe the HB 
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morphology of the majority of metal-rich GCs and thus they are of little use in the search for 
the second parameter. The solid error bars shown in Figure [11] represent the measurement 
errors listed in Table [TJ 

Metallicity uncertainty plays an important role in this diagram and thus the error bars 
in Figure [TT1 are only a lower limit to the total uncertainty in A(V — I). Unfortunately, it 
is not possible to provide realistic metallicity errors for all of the GCs in the sample. For 
the large deviation shown for —1.5 < [M/H] < —1 in Figure [TT] to be completely erased, 
metallicities of the GCs in this region that stand out would have to be systematically under- 
estimated by 0.25-0.5 dex. Such a conspiracy seems highly unlikely. However, it is possible 
to approximate the effect of metallicity errors on A(V — I) in the following manner. First, 
the fitting function described by equations [3] and @] was applied to the inner halo GCs with 
metallicities on the KI03 scale. Next, the KI03-based fitting function was used to calculate 
a A(V — I) but using the metallicity from Table [2] or the ZW84 scale, thereby introducing 
a metallicity error. Finally, the difference between the erroneous A(V — I) value obtained 
in this manner and the value obtained from the appropriate fitting function provides an 
estimate of how a metallicity error propagates into a A(V — I) error. 

Under the hypothesis that the KI03 [Fe/H] scale is the 'correct' one, this procedure 
gives an indication of how a metallicity error (a vertical displacement in Figure [TO]) in- 
fluences the significance of the distance between a GC and the inner halo fitting function 
(the horizontal distance between a point and the fit line in Figure [TO]) which we refer to 
as A (Fit — Measured). The dashed error bars in Figure [TT] are the quadrature sum of the 
A(V — I) measurement errors from Table [TJ and the metallicity errors estimated as described 
in the preceding paragraph. The metallicity effect is small except in the vicinity of the tran- 
sition from red to blue that occurs around [M/H] ~ — 1 (see Figure [TO]) where a small change 
in metallicity corresponds to a large change in A(V — I)- It should be stressed that this 
approach is not a rigorous treatment of the influence of metallicity uncertainty on A(V — I) 
or A (Fit — Measured) because it assumes the difference between one scale and another is 
the error. Nevertheless, it serves to illustrate that the existence of several GCs with —1.5 < 
[M/H] < —1 and large, positive values of A (Fit — Measured) cannot readily be attributed 
to metallicity errors alone. 

Figures [T2l through [T71 show how the A(V — I) difference depends on Rgc, My, Rh, 
R t , po, and age for the outer halo GCs. In order to quantify the information provided in 
the plots, Spearman's rank correlation coefficient was calculated for the data presented in 
Figures [T21 through [T71 and included in the figures. The Spearman rank correlation coefficient 
(Spearman's p) measures the degree of correlation between two variables but makes no 
assumption about the functional form of their relationship other than monotonicity. As 
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such, it is a very general measure of correlation; a perfect correlation has a value of p = +1, 
a perfect anti-correlation has p = — 1, and complete lack of correlation has p = 0. The left 
panels show the 6 most distant GCs, which were not part of the ACS Survey, as triangles. 
Since these 6 GCs are not part of the homogeneous ACS data set, and were not considered 
in Paper VII, we have calculated the Spearman correlation coefficients with and without 
them in the left panels. Apart from NGC 2419, the other 5 most distant GCs tend to 
strengthen the correlation coefficients because they are more distant, intrinsically fainter, 
more extended, and younger. 

A(V — I) error bars are included in Figures [10] through [T7J and age errors in Figure [TTj 
but uncertainties in the other quantities are neither readily available nor simple to estimate 
given the wide range of distances and reddenings present in the sample. It is, however, 
important to consider that the quantities in Figures [T21 through [TBI were converted to physical 
units after adopting a particular distance modulus. As distance increases, a given error in 
the distance corresponds to a larger uncertainty in the physical quantity. 

While Rqc, My, Rh, Rt> and p reveal no obvious visible trends, Figure [T7I shows a trend 
with age. The Spearman coefficients support this conclusion: My, R h , R t , and p show no 
strong correlations. Rgc shows a mild correlation but that is mostly likely due to a large 
anti-correlation between age and Rgc- Age shows a significant anti-correlation and the age 
trend works as anticipated: the HBs grow redder with decreasing age and so the distances 
between the measured A(V" — I) values and the inner halo trend increase as age decreases. 
The results presented in this section demonstrate that, of all the parameters considered here, 
age has the most significant correlation with A(Fit — Measured). 

Figure [T71 and the accompanying slopes and coefficients do not include Pal 12 or Terzan 7, 
the two youngest GCs in the sample. Both GCs are metal rich and so their HB morphologies 
are almost entirely determined by their metallicity alone. The presence or absence of Pal 12 
and Ter 7 does not significantly alter the correlation coefficients (Sp < 0.05) for any of the 
other quantities. Figure [17] also includes a linear fit to each data set that describes the rate 
at which HB morphology varies with age. The ages and metallicities listed in Table [2] suggest 
that HB morphology changes over a smaller age interval than those from Paper VII, but only 
by about ~0.5 Gyr. The larger errors that include the metallicity effect were not used in the 
least squares fits to the age- A (Fit — Measured) relations shown in Figure [T71 If they had 
been, the slope of the relation derived using the Paper VII ages (shown in the right panel) 
would be reduced to —0.21 ± 0.05 while the slope derived from the ages given in Table [2] 
would not change significantly. 
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5.2. A possible third parameter 



In §5.11 it was shown that when the inner halo relationship between metallicity and 
A(V— I) is subtracted from the outer halo GCs, a correlation between age and A (Fit — Measured) 
appears. The inner halo GCs were chosen to define the trend because they exhibit a tight 



relation with little scatter in the HB morphology-metallicity diagram (ISearle &: Zinnlll978 



see also Figure [TOl in this paper). Unfortunately, attempts to use this technique a second 
time-to subtract off the age trend- were unsuccessful because the remaining residuals are 
generally smaller than the measurement uncertainties in A(V — /). Thus, to identify a po- 
tential third parameter, it is necessary to remove the effect of metallicity and age to the 
fullest extent possible. The most metal poor GCs ([M/H] < —1.5) in the sample are an ideal 
choice because they have a small range in ages (2]) and only a weak metallicity dependence 
on HB morphology (see Figure [TO]) . 

Comparisons of the metal poor GCs' A(Fit — Measured) values with M v , R h , R t , and 
p were performed. Of these, central density (Log p ) produced the most significant trend. 
Figure [T8l demonstrates how central density relates to A (Fit — Measured) among the metal 
poor GCs using the two fitting functions employed in §5.11 and the same symbols used in 
Figures [TOl through [T71 There is an evident trend with the highest central density GCs having 
the bluest HB morphologies and hence the smallest A (Fit — Measured) values. Assuming 
that the contributions of age and metallicity are minimal in Figure dSJ the influence that 
central density has on HB morphology in the metal poor GCs is evident and its magnitude 
is ~ 0.2 in A(Fit — Measured). Although it is difficult to disentangle the effects of age and 
central density in GCs with intermediate metallicities (—1.5 < [M/H] < —1), where the 
second parameter is most pronounced, these GCs span essentially the same range of central 
densities as the metal poor GCs. Thus central density should only account for a small portion 
of the HB morphology variation seen in the intermediate metallicity group. The infl uence of 



central density o n HB m orph ology has alre ady been demonstrated by lFusi Pecci et al.l ( 119931 ). 
Buonanno et al.l (119971 ). and ISmithl (120041 ) . However, its effect is not as pronounced as that 
of age (compare Figures [T71 and [Tgj) . Therefore central density is the most likely candidate 
for the third parameter influencing HB morphology (as characterized by A(V — /)) that has 
been considered by this study. 



5.3. Discussion and comparisons with previous results 

It is worthwhile to demonstrate that the use of the inner halo GCs to define the HB 
morphology-metallicity relation does not include an unintended bias. In particular, since 
claims have been made in the preceding sections that, after metallicity, age and central 
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density are the two most significant factors influencing HB morphology [as represented by 
A(V — I)], it is important to demonstrate that the HB morphology-met allicity trend derived 
from the inner halo GCs does not include an implicit dependence on either of these quantities. 
Indeed, as Figure [TH] shows, the inner halo GCs exhibit some dispersion inside 8 kpc in both 
central density (left panel) and age (right panel) but there is no clear, systematic variation 
of either parameter as a function of Rgc inside 8 kpc (denoted by the dotted line in the 
figure). 

Figure UM demonstrates that the HB morphologies of the inner halo GCs, which have 
a fairly homogeneous age distribution, transition from blue to red as [M/H] rises through 

1. The metal rich GCs have red HBs (though NGC 6388 and NGC 6441 also have blue 

extensions) while the majority of metal poor GCs have blue HBs. It is at metallicities that 
are lower than, but within a few tenths of a dex of, [M/H]=— 1 where a GC is most likely to 
move a large distance in the A(V — I)— metallicity diagram due to a relatively small change 
in metallicity or another parameter. It is within this same narrow metallicity range that the 
AMR of the outer halo GCs branches off from that of the inner halo GCs (see Figure and 
so it is no coincidence that the classic second parameter effect is most pronounced in the 
outer halo. 

The analysis presented in §5.11 indicates that a GC with —1.5 <[M/H]< —1 requires 
2-2.5 Gyr to transition from a red HB to a blue one. This result is con sistent with a nu mber 
of studies mentioned in the introduction. We confirm the assertion by lLee et al.l (119941 ) that 
the presence of metal poor GCs with red HBs at large Rgc is primarily an age effect. How- 
ever, our result indicates that Lee et al.'s claim that the necessary age difference is greater 
than ~3 Gyr is slightly exaggerated^. That a GC with —1.5 < [M/H] < —1 transitions from 
a red to a blue HB in 2-2.5 Gyr poses a challenge to our understanding of mass loss during 
the red giant phase of evolution. RGB mass loss has be en essentially a free parameter in 
synthetic HB models since they were first constructed by iRoodl (119731 ). A number of mass 
loss rates that vary as a function of global stellar properties such as luminosity, mass, radius, 
or some combi nation of the three have been propo sed with va r ying degrees o f success, see 
Catelanl (120091 ) for a detailed discussion. Following iLed (Il99ll ). iDotterl (120081 ) showed that 
a simple relationship between global metallicity and average RGB mass loss, along with 
stochastic variations, can reproduce the general trend observed in old, roughly coeval Galac- 
tic GCs in the HB mor p hology- metallicity diagr a m. Recent prog r ess fr om ob servations by, 



e.g. lOriglia et all (120071 ). McDonald et all (120091 ). Meszaros et all (120091 ). and iDupree et al. 



7 It is important to keep in mind that there is a small, but not insignificant, number of GCs with younger 
ages, such as Pal 1 (Paper I), Pal 12, and Ter 7; these tend to be metal rich and so their HB morphologies 
are governed almost entirely by metallicity. 
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( 120091 ). should lead to a better understanding of RGB mass loss as a function of composition, 
evolutionary status, and pulsational properties. 



The recent, large-scale study of HB morphology by iRecio-Blanco et al.l (120061 ) suggests 
that, just as age can explain anomalously red HBs in some GCs, total GC mass is linked to 
the degree of abundance variation and the extent of the faint blue tail observed in other GCs. 
(Of course, there is no reason why both effects cannot be ope rating in the same GC.) R ecall 
that Figure [2] showed the HB morphology metric adopted by IRecio-Blanco et al.l (120061 ) was 
very sensitive to the extremes but lacked sensitivity through the middle of the distribution. 
In the absence of an ideal HB morphology metric, a large-scale study of HB morphology 
should choose a metric that is well-suited to the effect (s) of interest to that investigation. 



6. Conclusion 

HB morphologies characterized by A(V — I) of 66 Galactic GCs were examined to 
determine the sensitivity of HB morphology to a variety of different factors. Deep, homo- 
geneous photometry from the ACS Survey of Galactic GCs accounts for 60 of these while 
the remaining 6 GCs are the most distant Galactic GCs known. The complete sample is 
the largest examined to date and consists solely of high-quality HST photometry. It spans 
the full range of Rqc an d almost the entire range of metallicity present in the Galactic GC 
population. A(V — I) values and two independently measured sets of ages were joined with 
other quantities from the literature to assess the relative importance of these quantities on 
HB morphology. 

The data were split into two groups, roughly equal in number, based on Rgo The tight 
relationship between metallicity and HB morphology in the inner halo group (Rgc < 8 kpc) 
was characterized by a fitting function and this trend was subtracted off of the outer halo 
group. The difference between fit and measured A(V — /) was then compared to a variety of 
parameters, of which only age showed a significant correlation. The age correlation does not 
rely on the presence of the 6 most distant GCs in the analysis, though their presence does 
strengthen the result. Hence we conclude that, after metallicity, age has the most influence 
on A(y - 7). The age spread among the bulk of GCs in the Galactic halo was found to 
be 2-2.5 Gyr, though there are a few younger outliers such as Pal 12 and Ter 7. Further 
analysis, in which both metallicity and age were restricted, provided strong evidence that 
central luminosity density (po) is the third most influential parameter on A(V — I). 

We wish to express our gratitude to the anonymous referee for a thoughtful, thorough 
report that improved the accuracy and focus of the paper. We thank Ivan King for sharing 
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his differential reddening corrected CMDs in advance of publication and the following people 
for sharing their HST photometry: Peter Stetson (Pal 3, Pal 4, and Eridanus), Bill Harris 
(NGC 2419), and Eric Sandquist (NGC 2419). AD acknowledges support from a CITA 
National Fellowship and an NSERC grant to D. VandenBerg. Further support for this 
work (proposal number GO-10775) was provided by NASA through a grant from the Space 
Telescope Science Institute which is operated by the Association of Universities for Research 
in Astronomy, Incorporated, under NASA contract NAS5-26555. 
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F606W-F8 1 4W F606 W-F8 1 4W 

Fig. 1.— Portions of NGC 1261 and M 92 (NGC 6341) CMDs are shown. Thick lines 
indicate the regions that were selected to determine the median HB colors. Thin lines 
indicate measured quantities: vertical lines indicate the median colors of the HB and RGB; 
the horizontal line indicates the level of the HB. 
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Fig. 2. — Left: Comparison o f A(V — I) from this paper with (B— R)/(B+V+R) from 
Mackey fc van den Berghl (120051) for all 66 GCs in the sample. Right: Comparison of A(V— J) 
and Log T eff (HB) from iRecio-Blanco et al.l (120061 ) for the 30 GCs common to both studies. 
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Fig. 3. — Isochrone fits to the CMD of NGC 6362; the fit parameters are listed on the 
figure. The left panel shows the full CMD to indicate the agreement between the models 
and the data on the RGB and unevolved main sequence. The right panel focuses on the 
main sequence turnoff region to indicate how the age was determined. 
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Fig. 4. — Same as Figure [3] but showing the differential reddening corrected CMD of 
NGC 3201. 
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Fig. 5. — Same as Figure [3] but showing NGC 7099. 
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Fig. 6. — Absolute magnitude of the HB in F606W as a function of metallicity. The best 
fit line is shown as a solid line. The dashed lines indicate ±0.05 mag above and below the 
best fit line, which corresponds to the estimated uncertainty in the apparent magnitude of 
the HB as mentioned in §3.11 
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Fig. 7. — Comparison of the [Fe/H] scales from Table [2] and ZW84 with the KI03 scale for 
47 GCs. The horizontal and vertical scales are the same in both panels. The solid line in 
each panel indicates equality. 



-37- 




Fig. 8. — The normalized difference in age between this study and Paper VII as 
a function of metallicity. Age difference is presented in the sense that AAge = 
Age(ThisPaper)— Age(PaperVII) assuming, for consistency, ages from the latter based on 
the isochrones from Paper II and the ZW84 metallicity scale and normalized by dividing by 
the quadrature sum of the age uncertainties from both sources. 
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Fig. 9. — Top: Age-metallicity relations from the present study and Paper VII for the 55 
GCs in common. The open circles are the outer halo GCs. Bottom: Age-metallicity relations 
from the present study and IVandenBergl (120001 ) for the 20 GCs in common. 
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Fig. 10. — The metallicity-HB morphology diagram showing the inner (circles) and outer 
(squares) halo GCs from the ACS Survey. The six most distant GCs are shown as triangles. 
Error bars are from §3.2( fitting functions are shown as solid lines. The left panel shows the 
metallicity scale from Table [2J the right panels shows the ZW84 metallicity scale from Paper 
VII. 
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Fig. 11. — Difference between fitted and measured A(V — I) values as a function of [M/H]. 
The symbols and layout are the same as in Figure [101 The solid error bars are from mea- 
surement uncertainty in A(V — J) alone; the dashed error bars add the effect of metallicity 
uncertainty. 
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Fig. 12. — The difference between fit and measured A(V — I) values as a function of Rgc- 
Symbols and layout are as in Figure [TTJ The circles are M 5 (near the dotted line) and 
NGC 6584 (further away). Included on the figure is the value of the Spearman rank coef- 
ficient, see text for details. The phrase 'no OH' means that the 6 most distant GCs were 
omitted from the Spearman coefficient calculation. 
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Fig. 13. — Similar to Figure fT2l but as a function of My. 
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Fig. 14. — Similar to Figure [T2l but as a function of Rh- 
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Fig. 15. — Similar to Figure [T2l but as a function of R t . 
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Fig. 16. — Similar to Figure H2l but as a function of Log p . 
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Fig. 17. — Similar to Figure [12] but as a function of age. The left panel shows the ages 
reported in Table [2] while the right panel shows the D07/ZW84 ages from Paper VII. The 
solid line is a fit to the ACS Survey data alone; the dashed line is a fit that includes the six 
outer halo GCs. The two youngest GCs, Pal 12 and Ter 7, have been excluded from this 
figure and calculation of the coefficients, see text for discussion. 
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Fig. 18. — Central density (Log po) vs. A(Fit — Measured) showing only the metal poor 
GCs ([M/H] < —1.5). As with Figures [TT1 through \T7\ the left and right panels show results 
based on the two fitting functions determined in §5.11 Symbols are the same as in Figure [TUl 
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Fig. 19. — The dependence of central density (left panel) and age (right panel) on Rgc- The 
boundary assigned between the inner and outer halo is indicated by the dotted line at 8 kpc. 
The age plot shows results from both §4.21 (circles) and Paper VII (diamonds). Rgc an d Po 
data are from the Harris catalog. 
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Table 1. Horizontal Branch Parameters 







{V - J) HB 


(V - J) HB (V - /)rgb 




ID 


F606W(HB) 


Median 


MAD Median 


A(V-J) 



Arp 2 


18 


.05 


0.316 


± 


0.017 


0. 


.200 


1 


.056 


± 





.018 


0, 


.740 


± 





.025 


Lynga7 


17 


.10 


1.863 


± 


0.003 


0. 


.034 


1 


.974 


± 





.015 


0, 


111 


± 





.015 


NGC 104 


13 


.90 


0.916 


± 


0.002 





,023 


1 


.069 


± 





.002 





,153 


± 





.003 


NGC 288 


15 


.35 


-0.018 


± 


0.023 


0. 


.127 


1 


.004 


± 





.011 


1 


.022 


± 





.025 


NGC 362 


15 


.33 


0.786 


± 


0.003 


0, 


.083 





.981 


± 





.001 


0. 


,195 


± 





.003 


NGC 1261 


16 


.60 


0.762 


± 


0.004 


0. 


,129 





.963 


± 





.003 


0. 


,201 


± 





.005 


NGC 1851 


16 


.00 


0.764 


± 


0.011 


0. 


.272 





.998 


± 





.002 


0, 


,234 


± 





.011 


NGC 2298 


16 


.00 


0.332 


± 


0.018 


0. 


,123 


1 


.203 


± 





.007 





,871 


± 





.019 


NGC 2808 


16 


.10 


0.235 


± 


0.025 





,439 


1 


.201 


± 





.002 





,966 


± 





.025 


NGC 3201 


14 


.60 


0.767 


± 


0.108 


0. 


,270 


1 


.290 


± 





.011 


0, 


,523 


± 





.108 


NGC 4147 


16 


.85 


0.103 


± 


0.031 


0. 


,136 





.941 


± 





.011 





.838 


± 





.033 


NGC 4590 


15 


.65 


0.291 


± 


0.011 


0. 


,147 





.991 


± 





.003 


0, 


.700 


± 





.011 


NGC 4833 


15 


.50 


0.451 


± 


0.029 





,154 


1 


.351 


± 





.003 





.900 


± 





.029 


NGC 5024 


16 


.80 


0.090 


± 


0.009 


0. 


,088 





.941 


± 





.003 


0, 


,851 


± 





.010 


NGC 5053 


16 


.60 


0.145 


± 


0.041 


0. 


,158 





.927 


± 





.003 


0, 


.782 


± 





.041 


NGC 5272 


15 


.55 


0.207 


± 


0.014 


0. 


,208 





.943 


± 





.002 





,736 


± 





.014 


NGC 5286 


16 


.40 


0.372 


± 


0.011 


0. 


,182 


1 


.245 


± 





.003 


0, 


.873 


± 





.011 


NGC 5466 


16 


.45 


0.208 


± 


0.024 


0. 


,124 





.929 


± 





.006 


0, 


,721 


± 





.025 


NGC 5904 


14 


.90 


0.120 


± 


0.024 


0. 


,230 





.994 


± 





.001 





,874 


± 





.024 


NGC 5927 


16 


.35 


1.490 


± 


0.002 


0. 


,025 


1 


.602 


± 





.006 


0, 


.112 


± 





.007 


NGC 5986 


16 


.40 


0.357 


-|- 


0.011 


0. 


.136 


1 


,311 


± 





.004 





954 


± 


o 


.012 


NGC 6093 


16 


.20 


0.275 


± 


0.008 


0. 


,153 


1 


.198 


± 





.003 


0, 


.923 


± 





.008 


NGC 6101 


16 


.50 


0.203 


± 


0.007 


0. 


,091 


1 


.061 


± 





.005 


0, 


.858 


± 





.009 


NGC 6121 


13 


.17 


0.987 


± 


0.106 


0. 


.252 


1 


.528 


± 





.010 


0, 


.541 


± 





.106 


NGC 6144 


16 


.20 


0.612 


± 


0.009 


0. 


,100 


1 


.491 


± 





.007 


0, 


.879 


± 





.011 


NGC 6171 


15 


.40 


1.296 


± 


0.013 





.104 


1 


.499 


± 





.009 


0, 


,204 


± 





.015 


NGC 6205 


14 


.85 


-0.091 


± 


0.010 


0, 


.122 





.956 


± 





.002 


1 


.047 


± 





.010 


NGC 6218 


14 


.60 


0.227 


± 


0.009 


0. 


,099 


1 


.213 


± 





.004 


0, 


.986 


± 





.009 


NGC 6254 


14 


.80 


0.269 


± 


0.015 


0. 


.113 


1 


.290 


± 





.004 


1 


,022 


± 





.016 


NGC 6304 


15 


.95 


1.556 


± 


0.002 


0, 


.020 


1 


.662 


± 





.004 


0, 


.105 


± 





.004 


NGC 6341 


15 


.05 


0.036 


± 


0.013 


0. 


.093 





.946 


± 





.004 





,910 


± 





.013 


NGC 6352 


14 


.95 


1.249 


± 


0.004 


0. 


,025 


1 


.372 


± 





.013 


0, 


.123 


± 





.013 


NGC 6362 


15 


.18 


0.817 


± 


0.009 


0. 


.193 


1 


.064 


± 





.007 


0, 


.247 


± 





.012 


NGC 6366 


15 


.25 


1.817 


± 


0.010 





,030 


1 


.957 


± 





.011 





.141 


± 





.014 


NGC 6388 


17 


.00 


1.344 


± 


0.002 


0. 


,172 


1 


.470 


± 





.002 


0, 


.126 


± 





.003 


NGC 6397 


12 


.90 


0.217 


± 


0.007 





,066 


1 


.162 


± 





.010 





.944 


± 





.012 
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-( nnrinnori 










(V - J) HB 


(V - J) HB 


(V - /) RGB 




ID 


F606W(HB) 


Median 


MAD 


Median 


A(V-I) 


TV T /"I f < f* A A -t 

NGC 6441 


17.65 


1.501 ± 0.002 


0.154 


1.611 ± 0.002 


0.110 ± 0.003 


NGC 6496 


16.20 


1.224 ± 0.004 


0.023 


1.331 ± 0.007 


0.107 ± 0.008 


NGC 6535 


15.65 


0.554 ± 0.041 


0.127 


1.457 ± 0.011 


0.903 ± 0.042 


TV T /" A r* A 1 

NGC 6541 


15.15 


A 1 -1 A | AHA 

0.112 ± 0.010 


0.107 


1.076 ± 0.004 


0.964 ± 0.011 


NGC 6584 


16.40 


0.640 ± 0.061 


0.215 


1.048 ± 0.003 


0.408 ± 0.062 


NGC 6624 


15.85 


1.248 ± 0.003 


0.025 


1.383 ± 0.012 


0.135 ± 0.012 


NGC 6637 


15.75 


1.092 ± 0.002 


0.020 


1.230 ± 0.002 


0.138 ± 0.003 


NGC 6652 


15.77 


1.030 ± 0.004 


0.021 


1.171 ± 0.004 


0.141 ± 0.006 


NGC 6656 


14.15 


0.438 ± O.Oll 


0.141 


1.373 ± 0.004 


0.935 ± 0.012 


NGC 6681 


15.63 


0.074 ± 0.010 


0.112 


1.060 ± 0.002 


0.986 ± 0.010 


NGC 6717 


15.53 


0.268 ± 0.021 


0.128 


1.216 ± 0.003 


0.948 ± 0.021 


NGC 6723 


15.30 


0.704 ± 0.039 


0.301 


1.075 ± 0.003 


0.371 ± 0.039 


NGC 6752 


13.70 


-0.058 ± 0.009 


0.102 


1.024 ± 0.003 


1.082 ± 0.010 


NGC 6779 


16.15 


0.305 ± 0.007 


0.110 


1.220 ± 0.003 


0.914 ± 0.008 


NGC 6809 


14.35 


0.154 ± 0.021 


0.097 


1.060 ± 0.004 


0.906 ± 0.021 


NGC 6838 


14.21 


1.185 ± 0.009 


0.021 


1.291 ± 0.005 


0.106 ± 0.010 


NGC 6934 


16.78 


0.355 ± 0.046 


0.251 


1.066 ± 0.005 


0.711 ± 0.046 


NGC 6981 


16.73 


0.537 ± 0.122 


0.247 


1.008 ± 0.007 


0.471 ± 0.122 


NGC 7078 


15.75 


0.131 ± 0.028 


0.196 


1.007 ± 0.002 


0.876 ± 0.028 


NGC 7089 


15.95 


0.043 ± 0.009 


0.148 


0.993 ± 0.002 


0.951 ± 0.009 


NGC 7099 


15.12 


0.087 ± 0.005 


0.077 


0.960 ± 0.004 


0.872 ± 0.006 


Pal 12 


16.90 


0.908 ± 0.005 


0.008 


1.019 ± 0.016 


0.111 ± 0.017 


Terzan 7 


17.67 


1.028 ± 0.004 


0.012 


1.115 ± 0.005 


0.087 ± 0.007 


Tcrzan 8 


17.90 


0.230 ± 0.011 


0.111 


1.095 ± 0.011 


0.865 ± 0.015 


Additional GCs 


AM-1 


20.92 


0.799 ± 0.008 


0.021 


0.991 ± 0.005 


0.192 ± 0.009 


Eridanus 


20.23 


0.829 ± 0.015 


0.017 


1.009 ± 0.033 


0.180 ± 0.036 


NGC 2419 


20.35 


0.198 ± 0.008 


0.155 


1.017 ± 0.002 


0.819 ± 0.008 


Pal 3 


20.40 


0.688 ± 0.021 


0.087 


1.002 ± 0.019 


0.314 ± 0.028 


Pal 4 


20.65 


0.816 ± 0.008 


0.018 


0.999 ± 0.006 


0.183 ± 0.010 


Pal 14 


20.00 


0.828 ± 0.021 


0.035 


1.056 ± 0.027 


0.228 ± 0.034 



Note. — MAD = Mean Absolute Deviation 
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Table 2. Results from Isochrone Fitting 



Name [Fe/H] [a/Fe] DM F814W E(6-8) Age (Gyr) 



Arp 2 


-1 


.80 





.2 


17 


.55 


o 


.113 


13 


.00 


± 


0.75 




-0 


.60 


o 


,2 


15 


.80 





.713 


12. 


.50 


± 


1.00 


NGC 104 


-0 


.70 





.2 


13 


.30 


o 


.023 


12. 


.75 


± 


0.50 


NGC 288 


-1. 


.40 


o 


,4 


14 


.85 


o 


.013 


12. 


.50 


± 


0.50 


NGC 362 

X « V_ < V_.. KJ \J £-i 


-1. 


.30 





.2 


14 


.80 


o 


.023 


11. 


.50 


± 


0.50 


NGC 1261 


-1. 


.35 


0. 


,2 


16 


.10 


o 


.013 


11. 


.50 


± 


0.50 


NGC 2298* 


_1 


.90 


o 


.2 


15 


.43 


o 


.237 


13. 


,00 


± 


1.00 


NGC 3201* 


_1 


.50 


o 


.2 


13 


.90 


o 


.268 


12 


.00 


± 


0.75 


NGC 4147 


-1 


.70 


o 


,2 


16 


.48 


o 


.018 


12. 


.75 


± 


0.75 


NGC 4590 


-2 


.30 





.2 


15 


.30 


o 


.056 


13 


.00 


± 


1.00 


NGC 4833* 


_2. 


.30 


o 


,4 


14 


.84 


o 


.353 


13 


.00 


± 


1.25 


NGC 5024 


_2 


.00 





.2 


16 


.43 


o 


.023 


13 


.25 


± 


0.50 


NGC 5053 


_2. 


.40 


o 


,2 


16 


.32 


o 


.021 


13 


.50 


± 


0.75 


NGC 5272 


-1. 


.60 


o 


,2 


15 


.08 


o 


.018 


12 


.50 


± 


0.50 


NGC 5286* 


_1 


.70 


o 


.0 


15 


.75 


o 


.263 


13 


.00 


± 


1.00 


NGC 5466 


_2. 


.10 


o 


,2 


16 


.12 


o 


.023 


13. 


.00 


± 


0.75 


NGC 5904 


_1 


.30 





,2 


14 


.38 


o 


.033 


12. 


.25 


± 


0.75 


NGC 5927* 


-0. 


.50 


o 


,2 


15 


.30 


o 


.393 


12 


.25 


± 


0.75 


NGC 5986* 


_1 


.60 





,2 


15 


,73 





.295 


13 


.25 


± 


1.00 


NGC 6093* 


-1. 


.70 


o 


,2 


15. 


.55 


o 


.213 


13 


.50 


± 


1.00 


NGC 6101 


_1 


.80 


o 


,2 


16 


.03 





.113 


13 


.00 


± 


1.00 


NGC 6121* 


_1 


.20 


o 


,4 


12 


.20 


o 


.423 


12. 


.50 


± 


0.50 




i 

- 1 . 




n 
u 


o 

.z 


10 


/in 


u 


A AR 




.ou 


it 


1 nn 


NGC 6171* 


-1. 


.00 





.2 


14. 


.45 





.418 


12. 


.75 


± 


0.75 


NGC 6205 


-1. 


.60 





.2 


14. 


.47 





.019 


13 


.00 


± 


0.50 


NGC 6218 


-1 


.30 





,2 


13 


.90 





.191 


13 


.25 


± 


0.75 


NGC 6254* 


-1. 


.55 





.2 


14 


.15 


0. 


.261 


13 


.00 


± 


1.25 


NGC 6304* 


-0 


.50 





.2 


14. 


.87 





.473 


12. 


.75 


± 


0.75 


NGC 6341 


-2. 


.40 





.2 


14 


.80 





.031 


13 


.25 


± 


1.00 


NGC 6352* 


-0 


.80 





.4 


14 


.10 





.253 


13 


.00 


± 


0.50 


NGC 6362 


-1 


.10 





.4 


14 


.55 





.070 


12 


.50 


± 


0.50 


NGC 6366 


-0. 


.70 





.2 


14. 


.00 





.718 


12. 


.00 


± 


0.75 


NGC 6397 


-2. 


.10 





.2 


12 


.40 





.183 


13. 


.50 


± 


0.50 


NGC 6496 


-0 


.50 





.2 


15 


.35 





.213 


12 


.00 


± 


0.75 


NGC 6535 


-1. 


.90 





.2 


14 


.85 





.443 


13 


.25 


± 


1.00 


NGC 6541* 


-1 


.90 





.4 


14 


.68 





.118 


13 


.25 


± 


1.00 


NGC 6584 


-1. 


.40 





.2 


15 


.85 





.078 


12. 


.25 


± 


0.75 
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Table 2 — Continued 



Name 


[Fc/H] 


[a/Fe] 


DMpsuw 


E(6-8) 


Age (Gyr) 




-u.ou 


n n 
u.u 


10. uo 


n 9^1 


i o nn -4- n 7^ 


MPP 
INkjrL/ DDo / 


-U. i u 


n 9 

U.Z 


10. uo 


U.lOo 


1 o en _|_ n 7^ 


npp aaKO 
oooz 


-U. f 


9 

u.z 


10. uo 


U.llO 


i a ore _i_ n ^0 


NGC 6681 


-1.50 


0.2 


15.05 


0.098 


13.00 ± 0.75 


NGC 6717* 


-1.10 


0.2 


14.78 


0.203 


13.00 ± 0.75 


NGC 6723 


-1.00 


0.2 


14.67 


0.073 


12.75 ± 0.50 


NGC 6752 


-1.50 


0.2 


13.26 


0.053 


12.50 ± 0.75 


NGC 6779* 


-2.20 


0.2 


15.65 


0.248 


13.50 ± 1.00 


NGC 6809 


-1.80 


0.2 


13.88 


0.113 


13.50 ± 1.00 


NGC 6838* 


-0.70 


0.2 


13.40 


0.223 


12.50 ± 0.75 


NGC 6934 


-1.55 


0.2 


16.23 


0.108 


12.00 ± 0.75 


NGC 6981 


-1.50 


0.2 


16.20 


0.048 


12.75 ± 0.75 


NGC 7078 


-2.40 


0.2 


15.40 


0.083 


13.25 ± 1.00 


NGC 7089 


-1.60 


0.2 


15.48 


0.048 


12.50 ± 0.75 


NGC 7099 


-2.40 


0.2 


14.82 


0.053 


13.25 ± 1.00 


Pal 12 


-0.80 


0.0 


16.40 


0.033 


9.50 ± 0.75 


Terzan 7 


-0.60 


0.0 


17.15 


0.073 


8.00 ± 0.75 


Terzan 8 


-2.40 


0.4 


17.50 


0.133 


13.50 ± 0.50 


Additional GCs [DM V and E(V - /)] 


AM-1 


-1.50 


0.2 


20.41 


0.020 


11.10 ± 0.50 


Eridanus 


-1.30 


0.2 


19.78 


0.060 


10.50 ± 0.50 


NGC 2419 


-2.00 


0.2 


20.05 


0.062 


13.00 ± 1.00 


Pal 3 


-1.50 


0.2 


19.89 


0.060 


11.30 ± 0.50 


Pal 4 


-1.30 


0.2 


20.14 


0.055 


10.90 ± 0.50 


Pal 14 


-1.50 


0.2 


19.51 


0.045 


10.50 ± 0.50 



Note. — An asterisk (*) after the name indicates that the differential 
reddening corrected CMD was used in the isochrone analysis. E(6 — 8) 
= £(^606^-^814^). 
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Table 3. Fitting Function Coefficients 



Name 


This Paper 


Paper VII 


a 


0.947 ±0.016 


0.946 ±0.015 


ai 


0.809 ±0.017 


0.809 ±0.017 


a 2 


0.900 ±0.021 


1.012 ±0.036 


a 3 


0.022 ±0.011 


0.029 ±0.019 


bo 


-0.098 ±0.023 


-0.091 ±0.020 




-0.244 ±0.050 


-0.253 ±0.051 


b 2 


-0.105 ±0.020 


-0.127 ±0.025 



